Metallicity Gradients of Thick Disk Dwarf Stars 



Kenneth Carrell, Yuqin Chen and Gang Zhao 
Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy oj Sciences, 

Beijing 100012, China 
carrellQnao . cas . cn 



ABSTRACT 

We examine the metallicity distribution of the Galactic thick disk using F, G, and K dwarf 
stars selected from the Sloan Digital Sky Survey, Data Release 8. Using the large sample of 
dwarf stars with proper motions and spectroscopically determined stellar parameters, metallicity 
gradients in the radial direction for various heights above the Galactic plane and in the vertical 
direction for various radial distances from the Galaxy center have been found. In particular, 
we find a vertical metallicity gradient of -0.113 ± 0.010 (-0.125 ± 0.008) dex kpc^^ using an 
isochrone (photometric) distance determination in the range 1 < |Z| < 3 kpc, which is the 
vertical height range most consistent with the thick disk of our Galaxy. In the radial direction, 
we find metallicity gradients between +0.02 and +0.03 dex kpc~^ for bins in the vertical direction 
between 1 < |Z| < 3 kpc. Both of these results agree with similar values determined from other 
populations of stars, but this is the first time a radial metallicity gradient for the thick disk 
has been found at these vertical heights. We are also able to separate thin and thick disk stars 
based on kinematic and spatial probabilities in the vertical height range where there is significant 
overlap of these two populations. This should aid further studies of the metallicity gradients 
of the disk for vertical heights lower than those studied here but above the solar neighborhood. 
Metallicity gradients in the thin and thick disks are important probes into possible formation 
scenarios for our Galaxy and a consistent picture is beginning to emerge from results using large 
spectroscopic surveys, such as the ones presented here. 



Subject headings: galaxy: abundances — galaxy: disk 



1. Introduction 

The thick disk of our Galaxy has been exten- 
sively studied (especially in the solar neighbor- 
hood) since it was discovered bv lGilmore fc Reid 
(Il983h . Photometric and spectroscopic studies 



dominant component for stars with vertical dis- 
tances of 1 < \Z\ < 3 kpc f rom the Galactic plane 



(lAUende Prieto et al.lliooi ) 



have shown that it is a dist i nct co mponent in a 
number of ways. I Juric et "aL (2008) find that the 
stellar number density distribution of dwarf stars 
in the Galactic disk found using Sloan Digital Sky 
Survey (SDSS) data and a photometric parallax 
method is well fit with a two-component exponen- 
tial. The thick disk component is also known to 
have both a rotational lag an d a larger velocity dis- 
persion than the thi n disk ( Chiba fc Beerd 2000l : 
Soubiran et al ] l2003l) giving it a distinct kinematic 
signature. It has further been shown to be the 



The metallicity distribution of the thick disk 
is still an unresolved issue, however. Using 
the metallicity an d kinematics of SDSS stars, 
Ivezic et al. ( 20081) are able to reproduce the re- 



sults of previous works, but also caution that there 
are problems with describing their data using a 
two-component exponential. Most importantly, 
they suggest that the metallicity and velocity dis- 
tributions can be used to describe a single disk 
component instead of two separate ones. This 
would suggest a complex nature to a single disk 
component and not two distinct disk components 
with different formations and evolutions. More re- 
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centlv. lBovv et aL ( 2012 ) also find that the Galac- 
tic disk can be represented by a single exponen- 
tial in the radial and vertical directions that is 
smoothly varying as a function of age (by proxy 
using abundances). 

This has naturally lead to several compet- 
ing theories as to the formation mechanism of 
the thick disk and how it relates to the thin 
disk. For example, the thick disk could be 
mainly compo sed of accre t ed st ars from previ- 
ous mergers (|Abadi et al.l l2003l) 
stars for med from a ccreted 
mergers ([Brook et al 



2004 ) 



composed of 
gas from previous 
or from the inter- 
nal mixi ng of stars and g as within the Galaxy 
ISchonrich fc BinnevI l2009l ). These are just a 
few of the various internal and external sources 
from which the thick disk of our Galaxy could 
arise. Obviously, a combination of some or all of 
these scenarios could also give rise to our current 
Galactic hierarchy. Certain well-defined proper- 
ties of the thick disk should be able to distinguish 
between these possible scenarios and provide a 
clearer picture of the formation and evolution of 
our Galactic disk. 

In the radial direction, metallicity gradients 
vary from ~0 dex kpc~^ (lAllende Prieto et al 



2006 ) to -h0.028 ± 0.036 dex kpc'^ (jNordstrom et al 



20041 ) for stars likely to be thick disk members. 
We note that the former uses spectroscopically 
derived stellar parameters of F- and G-type stars 
observed by SDSS in the range 1 < \Z\ < 3 kpc 
while the latter uses photometrically derived stel- 
lar parameters of nearby F- and G-type stars. Re- 
cent res ults from the Radial Ve locity Experiment 
(RAVE, ISteinmetz et al.' 20o3) fall within these 



two values (|Ruchti et al., ,20 lit ICo^kunoglu et al 
I2OI2I) . 

Vertical metallicity gradie nts of the thick disk 
show more discrepancy. lAllende Prieto et al 



(|2006[ ) found no variation of metallicity as a func- 
tion of vertica l height. Using SD SS Red Horizontal 
Branch stars, IChen et al.l (|201lf ) found a gradient 
as high as -0.22 ± 0.07 dex kpc~^ but this re- 
sult most likely suffered contamination from other 
Galactic sources. Again, the RAVE results fall 



somewhere in the middle (|Kordopatis et al.ll2011 



Ruchti et al.ll2011) as do the results of Katz et al 



(|201l[ ) who observed stars along two different lines 
of sight. 

Our aim in this work is to examine both ra- 



dial and vertical metallicity gradients of stars as- 
sociated with the thick disk of our Galaxy. We 
use the large sample of spectroscopically observed 
stars from the Sloan Extension for Galactic Under- 
standing and Exploration (SEGUE) made avail- 
able in SDSS Data Release 8 (DR8). Selecting 
dwarf stars in this analysis is important since it is 
the dominant component of the disk of our Galaxy. 
Furthermore, these stars have changed very little 
since their formation, making them an excellent 
probe of the condition of our Galaxy during the 
formation of its disk. We employ various tech- 
niques used by previous analyses in order to more 
easily compare results from different sets of data 
and using different selection criteria. The large 
sample of stars with reliable parameter determi- 
nations facilitates this comparison and allows us 
to give a more comprehensive view of the metal- 
licity distribution of stars up to several kiloparsecs 
above or below the Galactic plane. 

2. Data 



Spectroscopic data from SDSS DR8 (jAihara et al 



20111 ) was obtained from the CasJobs interface. 
The improvement in DR8 for high-metallicity stars 
is important for this work since we expect many 
stars to have abundances near or slightly below 
solar. We used the seguel-targetl flag to select F, 
G, and K dwarfs (integer values of -2147483136, 
-2147221504, and -2147450880, respectively). We 
require that the stars have clean photometry and 
reasonable spectroscopic parameters (wide cuts 
were applied to eliminate any undetermined or 
extreme values for the effective temperature, sur- 
face gravity, metallicity and radial velocity) . This 
leaves us with 66,481 out of the total 88,081 stars 
for these three categories. Further cuts were ap- 
plied to ensure a clean sample of dwarf stars with 
good quality observed values and are summarized 
in TableHl Our final sample is 43,417 F, G, and K 
dwarf stars (slightly less than half the full sample) 
that have well defined photometric and spectro- 
scopic values as well as reasonable proper motions. 
In Figure [1] we show the spatial distribution of our 
dwarf star sample using distances calculated as 
described below. 

The most stringent cut for our dwarf sample 
is on the surface gravity since we want a sam- 
ple free of possible turnoff or giant stars to aid 
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Table 1: Selection cuts for our data 



Category 



SDSS Criteria 



Stars 



Our Criteria 



F/G 



G dwarf 



K dwarf 



g < 20.2 
0.2 <{g-r) < 0. 



14.0 <r < 20.2 
0.48 <{g-r) < 0.55 



14.5 < r < 19.0 
0.55 <{g-r) < 0.75 



6939 



62784 



18358 



12.0 < g,r,i < 21.0 mag 
g, r, i error < 0.05 mag 
log 5 > 4.4 dex AND log 5 error < 0.25 dex 
Teff error < 250 K 
[Fe/H] error < 0.25 dex 
radial velocity error < 10 km 
proper motion errors < 6 mas yr~^ 



TOTALS 



88,081 



43,417 




Fig. 1. — Spatial distribution of F, G, and K dwarf 
stars selected from SDSS DR8. 



in the reliability of the distance determination. 
Therefore, selecting a cut on the surface gravity 
of 4.4 dex and requiring that the error on this 
value be less than 0.25 dex ensures that we will 
not have stars with surface gravities lower than 
~4.15 dex. Furthermore, as can be seen in Figure 
[2l there is a noticeable change in the distribution 
of the metallicity as a function of surface gravity 
near log 5 ~ 4.4 dex. The histogram in the upper 
portion of Figure[2]clearly shows that the metallic- 
ity distribution of stars with 4.1 < logg < 4.4 dex 
is wider and has a peak that is more metal poor 
than for star s with log g > AA dex. We note that 
(|2012l ) concentrate on the main se- 



Cheng et al. 



quence turi ioff stars that we are tr ying to elimi- 
nate. Also, ISchlesinger et al.l (|201l[ ) use G and K 
dwarfs in their analysis, but their surface gravity 
cut is log (7 > 4.1 dex. This means that their sam- 
ple overl aps in surface gravi ties with this paper 
and with ICheng et al.l (j2012( ). The much steeper 



Fig. 2. — Surface gravities and metallicities for the 
F, G, and K dwarf stars selected from SDSS DR8. 
The horizontal lines in the lower plot correspond 
to surface gravities of 4.1 and 4.4 dex. The upper 
plot shows histograms of the metallicity for two 
surface gravity selections. 



vertical metallicity gradient in ISchlesinger et al 



(j201l| ) could be due to this fact as well as the mod- 
ifications they apply to the metallicities provided 
by the SEGUE Stellar Parameter Pipeline and the 
weighting algorithm they use (both of which tend 
toward higher metallicity). 

Selection bias in SDSS spectroscopic samples 
can be extremely difficult to determine or even es- 
timate. Our sample comes from a relatively simple 
set of selection criteria so we expect the bias to be 
minimal. Essentially, targets were selected from a 
wide color range (0.2 < {g ~ r) < 0.75) with more 
weight given to G type stars. Since we are not in- 
terested in a specific population, but in the dwarf 
stars as a whole, the number of each spectral type 
is not a primary concern. What is most essential is 
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that the selection criteria should not influence the 
expected abundances of our sample stars, and a 
simple color cut should not introduce a significant 
bias in this regard. 

Distance determinations were made using two 
methods. The firs t method we used was adapted 
from Equation 1 of lKordopatis et al.l j 2011 ), which 



Selected to match Ivezic et al . (2008) 



. 4 < (g-r) < 0.5 
. 5 < (g-r) < . 6 



is based on the method found in Zwitter et al 



( 2010t ). This technique uses theoretical isochrones 
to determine the absolute magnitude of a star with 
a certain set of parameters and their errors. This 
is, in essence, a Bayesian approach with the likeli- 
hood determined by the distance from a point on 
the isochrone and the prior probabilit y determined 



by the mass distribution. We used iGirardi et al 
(20041) isochrones that were downloaded from the 



CMl^ web interface. The grid of isochrones has 
ages of 2-12 Gyr in 1 Gyr increments and is linear 
in Z (logarithmic in [Fe/H]) from Z = 0.0001 to 
0.03 in 0.0001 increments. The range in [Fe/H] of 
roughly -2.25 < [Fe/H] < 0.2 dex includes the ex- 
pected metallicities of the dwarf stars in our sam- 
ple and the reliable limits of DR8. The large age 
range is necessary if we want to accurately deter- 
mine distances for both the (younger) thin disk 
and (older) thick disk stars, but since our sample 
is composed of only dwarf stars, this age range 
has little effect because all the stars in our sample 
should be main sequence stars. 

Using the mass weighted method described 
above, the probability distribution for each star 
was determined using four observables: the spec- 
troscopically determined effective temperature, 
surface gravity, and metallicity and the photo- 
metrically determined {g — r)o color. Individual 
errors on each parameter were used instead of 
overall errors quoted for DR8 measurements to 
give a better indication of the intrinsic uncer- 
tainty of the distance determination for each star. 
Once the probability distribution was calculated, 
the most probable absolute magnitude for the star 
was found from a weighted average and the error 
on this value was determined from the weighted 
standard deviation. Using this absolute magni- 
tude and the measured apparent magnitude one 
can find the distance to the star and convert to 
the Galactic Cartesian system (X,Y,Z) using its 
position on the sky. As a further quality control 




-1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 

Spectroscopic Metallicity (dex) 

Fig. 3. — Fractional difference between dis- 
ta nces determined u sing the isochrone method 
of iKordopatis et~aL (12011 ) an d the photometric 
method of llvezic et al.l (|2008l ) as a function of 
metallicity for two color ranges. 



^ http: / / stev.oapd.inaf.it / cgi-bin / cmd 



of our stellar sample, we only use stars that have 
distance errors from the isochrone method of less 
than 5%. This eliminates less than 0.1% of the 
sample, but does help identify and remove out- 
liers and stars with inconsistent observations or 
parameter determinations. 

The other distance determinatio n we us ed came 
from Equation 1 of Ivezic et al. ( 20081) . This 
method determines the absolute r magnitude of 
a star using its {g — i) color and metallicity. We 
use this method as a comparison to the isochrone 
method to help identify systematic errors in the 
distance determination. In Figure [3] we plot the 
fractional difference between the distances deter- 
mined using the isochrone and photometric meth- 
ods as a function of metallicity for two color selec- 
tions. In this plot we hav e used the same sele ction 
criteria as that used by llvezic et al. ( 20081 ) and 
determine distances using photometric metallici- 
ties derived using the equations from their work. 
The range of metallicities i n the plot is limite d to 
the reliable range found bv llvezic et al. I (fiooil) for 
their photometric abundance determination. We 
find that the distan ces derived using the method of 
Ivezic et al.l (|2008l ) are sy stematically farth e r than 
the isochrone method of lKordopatis et aD (2011) 
by ~10% on the high metallicity end ([Fe/H] ^ - 
0.55 dex) and increase to '^20% at lower metal- 
licity ([Fe/H] ~ -1.15 dex). If, however, we use 
the spectroscopically determined [Fe/H] values in- 
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stead of the photometric ally determined on es in 



the distance estimation of llvezic et alJ (|2008l ). we 



find a nearly constant systematic difference of 
~10% for the entire metallicity range shown in 
Figure [31 

Using the distance, radial velocity, and proper 
motion of each star, the full space motions 
(U,V,W) were determ i ned u sing the method in 



Johnson fc SoderblomI ([1987|) with an updated 
galactic coordinate transformation. We used a 
modified version of the gaLuvw routine from the 
IDL Astronomy User's Librar^ and th e solar mo- 
tion values from lSchonrich et al. ( 201Clh to correct 
to the local standard of rest. The uncertainties 
in each of the velocity components are on average 
~20 km s~^ for our sample of stars. 

3. Analysis 

The cuts described in the previous section were 
applied to provide a sample of dwarf stars that 
could be used to isolate likely members of the disk 
component (s) of our Galaxy. As was described in 
Section [H past studies have used various methods 
to select and characterize the disk. We restrict our 
samples to 7.0 < R < 10.5 kpc since these radial 
distances provide us with the largest statistics and 
to avoid possible contamination from the bulge of 
the Galaxy. In the vertical direction we limit the 
sample to 0.5 < \Z\ < 3.0 kpc for determining ra- 
dial metallicity gradients since there are few stars 
below these heights and the statistics fall rapidly 
above them. For vertical metallicity gradients we 
fit in the range 1.0 < \Z\ < 3.0 kpc to be more con- 
sistent with previous works and for reasons that 
will be discussed in Section U) These restrictions 
in radial and vertical distances are shown in Figure 
[l]with solid lines. 

To determine the extent to which the thin disk 
contributes to our dwarf star sample we use the 
kinematic properties of the different components 
of the Galaxy as a discriminant. To do this, we 
used Equation 1 of iBensbv et al. ( 2003 ) and the 
characteristic velocity dispersions and asymmet- 
ric drift values given in their Table 1 to determine 
the likelihood for each star of belonging to either 
the thin disk, thick disk, or halo based on its kine- 
matics alone. The likelihood for each star of be- 



longing to a particular component of the Galaxy 
also depends on the local num ber density expe cted 
of that component. Since .Bensbv et al ] (l2003h use 
local densities for the solar neighborhood and our 
sample of stars is at higher v ertical heigh t s, we use 
the Galactic model given by I Juric et al. ( 2008 ) in 
column 2 of their Table 10. This is their best-fit 
model determined directly from the number den- 
sity distribution maps found in their work. Com- 
bining the kinematic and local number density 
likelihoods, for each individual star we determine 
the likelihood of it belonging to the thin disk, thick 
disk, or halo. Ratios of these likelihoods were then 
used to find stars that are more likely to be thin 
disk than thick disk (i.e. TD/D < 0.10) and more 
likely to be thick disk than thin disk or halo (i.e. 
TD/D > 10 and TD/H > 10). These limits of 
being ten times more likely is in agreement with 
other papers employing the same technique (e.g. 



Coskunoglu et al.ll2012l) . This helps separate stars 



^http: / /idlastro. gsfc.nasa.gov/ 



with kinematics most obviously matching either 
the thin or thick disks at the heights in which both 
populations overlap. We note that systematics in 
our distance determinations can affect the local 
number densities for each star and can also lead 
to incorrect velocity components, and therefore in- 
correct likelihoods for belonging to one Galactic 
component or another. This possible issue is mit- 
igated by the fact that we use strong cuts on the 
populations (i.e. being ten times more likely to 
be one than the other) and only use the likeli- 
hoods to separate our sample at the smallest ver- 
tical heights. 

In the lower panel of Figure H] we show that the 
vertical distributions for stars defined as likely thin 
disk and likely thick disk are very different. Thin 
disk stars are peaked at heights of slightly less than 
0.5 kpc and have no stars with heights larger than 
~0.75 kpc. Meanwhile, the thick disk distribution 
is much broader and has a peak slightly higher 
than 1 kpc. The upper plot in Figure |4] shows 
that the kinematically selected thick disk sample 
has a peak metallicity near -0.6 dex, while the 
thin disk sample is peaked at a higher metallicity 
(near -0.2 dex). Both the vertical height distribu- 
tion and the mean metallicities of the kinemati- 
cally selected components are consistent with pre- 
vious works. Of note is the fact that it is possible 
to separate the thin and thick disks in SDSS data 
using their kinematics and distances. Also, we see 
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0.5 1 1.5 2 2.5 3 3.5 

|Z| (kpo) 

Fig. 4. — Distribution in \Z\ (lower) and [Fc/H] 
(upper) for stars selected according to probabili- 
ties from spatial and kinematic information as de- 
scribed in the text. 

that for vertical heights of 0.5 < \Z\ < 1.0 kpc 
our dwarf star sample contains both thin and thick 
disk stars (~80% of the kinematically selected disk 
stars arc thick disk), while above \Z\ = 1 kpc we 
only have thick disk stars. 

Using the large sample of dwarf stars with reli- 
ably determined properties and a number of pre- 
vious values with which to compare, we can now 
provide a better view of the overall abundance pat- 
terns of dwarf stars in the Galactic thick disk in 
both the radial and vertical directions. 

4. Results 

When looking for metallicity gradients of the 
disk many factors arc involved that could affect 
the outcome. The population of stars, their lo- 
cations in the Galaxy, and their kinematics are 
all distinguishing characteristics that can (and do) 



give different results. In this paper we use exclu- 
sively dwarf stars that are mostly from the thick 
disk component of our Galaxy. We have a very 
large sample, however, with a wide range of dis- 
tances and kinematics. As such, we divide our 
sample based on these properties. Specifically, as 
discussed in the previous section, we find that for 
vertical heights oi \Z\ > 1 kpc our sample has no 
thin disk contamination and for 0.5 < \Z\ < 1.0 
kpc the contamination is approximately 20%. Fur- 
thermore, we group our stars by their positions in 
order to disentangle the affect of one metallicity 
gradient on the other. For example, if wc sim- 
ply took our entire sample and fit for a vertical 
metallicity gradient using all radial distances, the 
contribution of a radial metallicity gradient could 
have an affect on the vertical metallicity gradient. 
We choose bin sizes of 0.5 kpc in order to have a 
statistically significant sample of stars in each bin 
and minimize the affect of the other gradient. We 
also provide metallicity gradients without binning 
to see how these values agree with those previously 
published. 

It is important to note that while we do have 
stars in the full range of fitted values, not all 
the bins will have the same distribution of radial 
and/or vertical distances. This could affect our 
gradients. By averaging the points in each bin, 
however, and then determining the gradients from 
these averaged values we reduce the affect of these 
changes and any potential outliers in our sample. 

To estimate errors in the metallicity gradients, 
a bootstrapping method was employed where we 
randomly re-sampled half the number of stars in 
each group one thousand times and fit each one for 
a gradient. The errors quoted in the tables that 
follow are the standard deviations of the results 
from this method. To determine the affect of a 
systematic difference in the distance estimates we 
have included metallicity gradients using both the 
isochrone and photometric distance methods. 

4.1. Radial Gradients 

We first look for metallicity gradients in the 
radial direction for our sample stars. The pres- 
ence of a radial metallicity gradient in the Galac- 
tic plane has been observed using varioiis meth- 
ods (blue stars, open clusters, nebulae, etc). A 
radial metallicity gradient away from the Galactic 
plane is not as well established. In fact, depend- 
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ing on the population and he ight above the plane 



differe nt results are found. lAllende Prieto et al 
( 20061 ) find no evidence of a radial metallicity gra- 
dient for stars in the range 1 < |Z| < 3 kpc and 
while iNordstrom et alj (|2004[ ) find a positive gra- 
dient for their old age sample (>10 Gyr), their 
value is still consistent with no gradient within er- 
rors. Recent results from RAVE have found posi- 
tive gradients for likely thick disk stars selected in 
various ways and their results are consistent with 
a positive, nonzero gradient within errors. Radial 
metallicity gradients from various references can 
be found in Table El 

Since we have such a large sample, we are able 
to look not only for a radial metallicity gradient, 
but also see how the gradient changes as a func- 
tion of height above the Galactic plane. In Figure 
[5] we show metallicity as a function of the radial 
distance from the Galactic center for our entire 
sample of stars for different ranges in the vertical 
height. The resulting radial metallicity gradients 
can be found in Table [3] and a plot of the radial 
metallicity gradient as a function of average \Z\ 
is shown in Figure [6l For reference, horizontal 
lines have been drawn in Figure H] for the results 
of previous publications found in Table [21 and two 
overlapping points in \Z\ from lCheng et all (|2012[) 
have been included as well. For vertical heights of 
1 < \Z\ < 3 kpc we find that the radial metallicity 
gradients are relatively constant as a function of 
\Z\. In the vertical height bin where thin disk stars 
contaminate our sample (0.5 < \Z\ < 1.0 kpc) 
the radial metallicity gradient is slightly lower. In 
the highest vertical height bin (2.5 < \Z\ < 3.0 
kpc) the number of stars is an order of magnitude 
smaller than the bins with lower vertical heights 
and the difference between the results using the 
two different distance estimates is largest. Fur- 
thermore, from Figure [5] one can see that in this 
vertical height range the point for radial distances 
of 10.0 < R < 10.5 kpc using the isochrone method 
is noticeably above the fit and forces the linear 
fit to be steeper. Disregarding this point when 
fitting for the metallicity gradient gives a value 
(-1-0.023 dex kpc^^) that is more consistent with 
lower heights and the photometric distance esti- 
mate. We have included this value in Figure [S] as 
the open square. 

We could fit for a radial metallicity gradient 
for the entire range associated with the thick disk 
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Fig. 5. — Fits for the radial metallicity gradient 
of the dwarf star sample in different \Z\ ranges. 
Contour hues are shown for 95%, 90%, 68%, 50%, 
33%, and 10% of the peak density. 



(1 < |Z| < 3 kpc), however, as is evident from 
the offsets in the linear fits shown in Figure \5\ 
there is a strong vertical metallicity gradient in the 
sample that will affect an overall determination in 
the radial direction. 

Radial metallicity gradients determined using 
the photometric method distances are very simi- 
lar to those found with the isochrone method. In- 
deed, the values all agree with each other within 
errors. The fact that there is no large discrep- 
ancy between the results using the two methods 
means that systematic differences in the distances 
on the order we observe between the two distance 
estimates (up to ^20%) do not significantly affect 
our results. 

Our radial metallicity gradients lie in the mid- 
dle of the previously pu blished values and agree 
most closely with those of INordstrom et ah ( 2004 ) 
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Table 2: Radial metallicity gradients from the literature. 





gradient (dcx kpc ^) 


error 


notes 


Allende Pricto ct al. (2006) 


-0 




1< \Z\ <3 kpc 


Nordstrom ct al. (2004) 


+0.028 


0.036 


age >10 Gyr 




+0.016 


0.011 


F dwarfs, TD/D>10 


Coskunoelu et al. (2012~) 


+0.016 
+0.010 


0.012 
0.009 


F dwarfs, e.^ >0.1 
G dwarfs, TD/D>10 




+0.037 


0.016 


G dwarfs, Ct, >0.1 


Ruchti et al. (2011) 


+0.01 


0.04 


[Fe/H]<-1.2 


ChenE et al. (2012) 


+0.0028 
-0.013 


+U.I)UVi 
-0.0052 
+0.0093 
-0.0016 


1< \Z\ <1.5 kpc 
0.5< \Z\ <1.0 kpc 



Table 3: Radial metallicity gradients of our dwarf star sample. 





Isochrone Distance Method 


Photometric Distance Method 


\Z\ Height 


N 


\z\ 


R 


Gradient 


N 


\z\ 


R 


Gradient 


(kpc) 


stars 


(kpc) 


(kpc) 


(dcx kpc^^) 


stars 


(kpc) 


(kpc) 


(dex kpc^^) 


0.5-1.0 


16,892 


0.74 


8.46 


+0.015+0.005 


15,150 


0.75 


8.52 


+0.017+0.005 


1.0-1.5 


8,881 


1.22 


8.47 


+0.025+0.006 


9,424 


1.22 


8.50 


+0.019+0.006 


1.5-2.0 


4,314 


1.72 


8.53 


+0.025+0.007 


4,954 


1.72 


8.56 


+0.021+0.007 


2.0-2.5 


1,907 


2.22 


8.59 


+0.029+0.010 


2,519 


2.22 


8.56 


+0.027+0.008 


2.5-3.0 


730 


2.71 


8.68 


+0.041+0.016*^ 


1,136 


2.72 


8.64 


+0.022+0.014 



"Excluding the point with radial distances of 10.0 < R < 10.5 kpc gives a metallicity gradient of +0.023 dex kpc . 



and t he F dwarf samples of ICoskunoglu et al 



( 2012 ). Our results differ from those of lCheng et al 



20121 ) in the bins that overlap with their work. As 



is shown in Figure [U there is a noticeable change 
in the metallicity distribution for stars with sur- 
face gravities below ~4.4 dex in our sample, which 
could be the reason for the discrepancy between 
our work using high surface gravity dwarfs and 
their work using turnoff stars with lower surface 
gravities. It is also possible that there are differ- 
ent ratios of thin and thick disk stars between our 
sample and theirs, which would obviously give a 
different metallicity gradient if the two popula- 
tions themselves had gradients different from each 
other. 

4.2. Vertical Gradients 

In a similar fashion to what was done for the ra- 
dial gradients, we examine the vertical metallicity 
gradients of our sample. In Table |4] we list some 
previously published values using various popula- 
tions and sample selection methods. 



We bin in radial distance to fit for vertical 
metallicity gradients and restrict our fits to verti- 
cal heights of 1 < |Z| < 3 kpc for two reasons. The 
first is to be more consistent with previously pub- 
lished values. The second is to provide a clearer 
picture of the vertical metallicity gradient of only 
the thick disk. There is a significant overlap of 
the thin and thick disks for vertical heights of 
\Z\ < 1 kpc, which can be seen in the vertical dis- 
tribution of the kinematically separated samples 
shown in Figure |4l Therefore, if we extend our 
fits below this height for our full sample of stars 
we will not be probing a single population of our 
Galaxy, but the contribution from both the thin 
and thick disks. As we will show, the kinemati- 
cally selected thick disk with \Z\ < 1 kpc could 
be included in the fits without significantly alter- 
ing the results, but in order to be more consistent 
with previous results we use only vertical heights 
of 1 < |Z| < 3 kpc. 

Vertical metallicity gradients for our sample of 
stars are given in Table [5l Plots showing the verti- 
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Table 4: Vertical metallicity gradients from the literature. 





gradient (dcx kpc ^) 


error 


notes 


Allende Prieto et al. ('2006) 


-0 




1 < |Z| < 3 kpc 


Chen et al. C2011) 


-0.120 
-0.225 


0.01 
0.07 


Galaxy model and fits 
All points 


Katz et al. (2011) 


-0.068 


0.009 


2 lines of sight 


Ruchti ct al. (2011) 


-0.09 


0.05 


[Fe/H]<-1.2 


Kordopatis ct al. (2011) 


-0.14 


0.05 


1 < |Z| < 4 kpc 
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Fig. 6. — Radial metallicity gradients of our 
dwarf star sample found using distances from the 
isochrone (solid squares) and photometric (solid 
circles) methods for different heights above the 
Galactic plane. The horizontal lines and open tri- 
angles correspond to previously published values 
(see Table [21). 



cal metallicity gradients in different radial bins can 
be found in Figure [71 This plot shows the effect of 
selecting only stars with heights above \Z\ — 1 kpc. 
The points in the 0.5 < \Z\ < 1.0 kpc bin show 
the average metallicity of the full sample (open 
squares) and of the kinematically selected thick 
disk sample (open circles). Careful examination 
of these points shows that the full sample has a 
higher metallicity and is systematically above the 
fit value. This is most certainly due to the con- 
tribution of the more metal-rich thin disk in our 
sample at these heights and is the reason we have 
excluded this region from our metallicity gradient 
determinations. The kinematically selected thick 
disk sample matches well with the fit value. 

A plot of the vertical metallicity gradients as 
a function of radial distance from the Galactic 



center is found in Figure [3 It is obvious that 
the metallicity gradient for the largest radial dis- 
tances are the most discrepant between the two 
distance determinations. As was noted for the ra- 
dial metallicity gradients, these distances have the 
lowest number of stars making statistics a pos- 
sible issue in the metallicity gradient determina- 
tion. From Figure [3 one can see that including 
the 0.5 < \Z\ < 1.0 kpc bin in the fit for radial 
distances of 10.0 < R < 10.5 kpc found using 
the isochrone method would give a steeper slope. 
This would bring the result closer to the values at 
smaller Galactic radii as well as the result using 
photometric distances. As was stated above, how- 
ever, for consistency with previous results we have 
restricted our fits to vertical heights of 1 < \Z\ < 3 
kpc. 

There is some structure to the gradients as a 
function of radial distance but no obvious trends. 
In fact, a single value can describe our results. 
Because the radial metallicity gradient is much 
smaller (by an order of magnitude) than the ver- 
tical metallicity gradient, we can fit the data for 
stars with radial distances of 7.0 < R < 10.5 kpc 
and vertical heights of 1 < \Z\ < 3 kpc, which gives 
us values of -0.113 ± 0.010 and -0.125 ± 0.008 
dex kpc~^ for our full dwarf star sample using 
isochrone and photometric distances respectively. 
These values are in very good agreement with 
the shallower vertical metallicity gradient found 
bv IChen et all (1201 ll) as w ell as those g i ven i n 



Kordopatis et al.l (|2011l ) and lRuchti et all (1201 ih 
It is s lightly steeper than the result of iKatz et al 
(|201ll) . 



Vertical metallicity gradients found using the 
photometric distance determination are very simi- 
lar to those found using the isochrone distance de- 
termination. The values from each method agree 
within errors and there is no obvious systematic 
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Table 5: Vertical metallicity gradients of our dwarf star sample with heights of 1 < |Z| < 3 kpc. 





Isochrone Distance Method 


Photometric Distance Method 


R Distance 


N 


\z\ 


R 


Gradient 


N 


\z\ 


R 


Gradient 


(kpc) 


stars 


(kpc) 


(kpc) 


(dex kpc^^) 


stars 


(kpc) 


(kpc) 


(dex kpc"-^) 


7.0-7.5 


1,812 


1.58 


7.27 


-0.13±0.03 


2,042 


1.67 


7.26 


-0.12±0.03 


7.5-8.0 


3,125 


1.52 


7.76 


-0.09±0.03 


3,473 


1.57 


7.76 


-0.12±0.02 


8.0-8.5 


3,427 


1.48 


8.24 


-0.12±0.02 


3,827 


1.52 


8.25 


-0.13±0.02 


8.5-9.0 


2,913 


1.54 


8.74 


-0.16±0.02 


3,286 


1.57 


8.74 


-0.16±0.02 


9.0-9.5 


2,296 


1.55 


9.23 


-0.11±0.02 


2,674 


1.59 


9.23 


-0.13±0.02 


9.5-10.0 


1,423 


1.65 


9.73 


-0.11±0.02 


1,684 


1.67 


9.73 


-0.09±0.02 


10.0-10.5 


836 


1.67 


10.24 


-0.07±0.03 


1,047 


1.74 


10.24 


-0.13±0.03 


7.0-10.5 


15,832 


1.54 


8.51 


-0.113±0.010 


18,033 


1.59 


8.53 


-0.125±0.008 



difference between the two. The only discrepancy 
between the metallicity gradients determined us- 
ing the different distance estimates is for the bin 
with 10.0 < R < 10.5 kpc, which also has the 
lowest statistics. 

5. Discussion 

We have probed the metallicity gradients in 
both the radial and vertical directions using dwarf 
stars selected from SDSS. We kinematically sep- 
arate this sample into thin and thick disk com- 
ponents at vertical heights where the two popu- 
lations are known to overlap using probabilities 
defined and used in previous analyses. This is sig- 
nificant because this separation is most often ap- 
plied to stars in the solar neighborhood, but we 
have shown that it can be successfully used for 
SDSS data with vertical heights of up to 1 kpc by 
taking into account the local number density as a 
function of height. 

Stars associated with the Galactic thick disk 
have a positive radial metallicity gradient for 
heights of 1 < |Z| < 3 kpc. Our results seem 



to be at odds with the results of iCheng et al 



(|2012l ) in the two vertical height bins that overlap. 
However, the different populations and selection 
methods between our work and theirs make it dif- 
ficult to compare the two results directly. Our 
results indicate that a radial metallicity gradient 
determined for stars with \Z\ < 1.0 kpc will heav- 
ily depend on the fraction of thin and thick disk 
stars in the sample if the two populations have 
very distinct values and trends. 



When fitting the entire sample of stars with 
vertical distances most consistent with the thick 
disk, we find a radial metallicity gradient that is 
positive and not consistent with zero within er- 
rors. In particular, we find radial metallicity gra- 
dients of -1-0.025 ± 0.006 (+0.019 ± 0.006), -f 0.025 
± 0.007 (4-0.021 ± 0.007), +0.029 ± 0.010 (+0.027 
± 0.008), and +0.041 ± 0.016 (+0.022 ± 0.014) 
dex kpc~^ using isochrone (photometric) distances 
for stars in 0.5 kpc bins in vertical height between 
1 < \Z\ < 3 kpc and with radial distances of 
7.0 < R < 10.5 kpc. Restricting the fit for the 
highest vertical height bin to 7.0 < R < 10.0 kpc 
changes the result using isochrone distances from 
+0.041 to +0.023 dex kpc~\ which agrees bet- 
ter with lower heights and with the photometric 
distance determination. To our knowledge, this 
is the first time that a radial metallicity gradi- 
ent for the thick disk has been found at these 
vertical heights. The values we find are in very 
good agree ment wit h the old stellar sample value 
of Nor dstrom et al. 1 12004) and with both of the 



thick disk F dwarf samples of iCo^kunoglu et al 
(j2012r ). It al so agrees with in errors with the re- 
sults of iRuchti et al.i (|201l[) . Due to the size of 



our sample, our statistical error is smaller than 
these previous works. 

Fitting the full sample of stars with radial dis- 
tances of 7.0 < R < 10.5 kpc in the range 1 < 
|.Z| < 3 kpc, we find a vertical metallicity gradi- 
ent of -0.113 ± 0.010 (-0.125 ± 0.008) dex kpc"! 
using isochrone (photometric) distances. This 
value is in good a greement with the resu lts of 
(|201ll ). iKordopatis et al.l (|201l[ ). and 



Chen et al 
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Fig. 7. — Fits for the vertical metallicity gradi- 
ent of the dwarf star sample in different R ranges. 
Contour lines are shown for 95%, 90%, 68%, 50%, 
33%, and 10% of the peak density. 



Ruchti et all (|20m 

The results of Bow et al. ( 20121) (among oth- 
ers) suggests that the Galactic disk is a single 
component that changes smoothly as a function 
of both radial and vertical distances. We observe 
no obvious or significant trends in either the radial 
or vertical metallicity gradients for stars most con- 
sistent with the thick disk of our Galaxy. However, 
it is not clear if the kinematic selection, distances, 
and/or bin sizes we have used would be able to 
detect changes in the metallicity gradient due to 
such a scenario. Furthermore, as discussed in their 
work, dividing the disk into components using dif- 
ferent methods can give conflicting results, making 
it difficult to compare the results from geometric, 
kinematic, and abundance decompositions. 

In summary, the metallicity gradients we find in 
both the radial and vertical directions using the 



8 . 5 

R (kpc) 



Fig. 8. — Vertical metallicity gradients of our 
dwarf star sample found using distances from the 
isochrone (solid squares) and photometric (solid 
circles) methods for different radial distances from 
the Galactic center. The horizontal lines corre- 
spond to previously published values (see Table 
15). 

sample of dwarf stars selected from SDSS DR8 
show no large or obvious trends in the distance 
ranges we are able to probe. The vertical metal- 
licity gradient we have found is much steeper than 
the radial metallicity gradient, but we have found 
both to be nonzero. Future large sp ectroscopic 



surve y programs, such as LAMOST (jZhao et al 
l200(3l ). should provide greater statistics and more 
complete samples of stars in a wider range of dis- 
tances so that metallicity gradients can be probed 
in more detail and an even better understanding 
of our Galactic disk can be realized. 
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